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Abstract 

We present a simple mathematical model in which a time averaged pattern 
emerges out of spatio-temporal chaos as a result of the collective action of 
chaotic fluctuations. Our evolution equation possesses spatial translational 
symmetry under a periodic boundary condition. Thus the spatial inhomo- 
geneity of the statistical state arises through a spontaneous symmetry break- 
ing. The transition from a state of homogeneous spatio-temporal chaos to 
one exhibiting spatial order is explained by introducing a collective viscosity 

which relates the averaged pattern with a correlation of the fluctuations. 
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Macroscopic structures formed in turbulent flow, such as roll patterns seen in clouds 
and the red spot on Jupiter, are often observed in nature. These phenomena have been 
explained using a plausible argument that an equation describing mean flow is renormalized 
by turbulent fluctuations and displays an instability under certain conditions. Much progress 
has been made toward justifying this explanation [|T],0], but there has been little study of the 
origin of the mean flow. If a system possesses spatial inhomogeneity, it is easily understood 
how mean flow can arise, but this is not the case in which we are interested. In this paper, 
we consider the emergence of mean flow through a spontaneous symmetry breaking of the 
statistical state in turbulence. 

Several recent experiments seem to show the formation of a mean flow in turbulence. In 
thermo-convective systems, two different types of turbulence, soft and hard turbulence, can 
be observed, depending on the Rayleigh number, and it has been found that hard turbulence 
is characterized by the existence of mean flow . Because the experiments in question have 
been done on systems with small aspect ratios, effects of horizontal boundaries cannot be 
neglected, and the spatial translational symmetry is broken externally. However, the systems 
do possess discrete symmetry associated with a spatial parity. The mean flow then breaks 
the parity symmetry spontaneously, that is, the direction of the large scale circulation flow is 
determined by the initial conditions. We note that according to recent numerical simulations 
of the Boussinesq equation the horizontal boundaries do not play an essential role in the 
formation of mean flow. We thus expect that mean flow may occur even in the infinite limit of 
the aspect ratio. Actually, the emergence of mean flow in convective turbulence in a system 
with large aspect ratio (about 1000) has been found in the electro-hydrodynamic convection 
of nematics 0, where persistent roll patterns appear in the turbulent state. In this case, 
the mean flow can be thought to spontaneously break the spatial translational symmetry in 
turbulence. Additionally, it has been reported recently that chaotic wave patterns produced 
by the Faraday instability display spatially ordered time averages 0. Similar phenomena 
have also been discovered in Taylor-Coutte flow and rotating convection [Q. 

These experimental results suggest that mean flow and persistent patterns can appear 
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spontaneously in various chaotic systems. We propose to call these 'collective patterns in 
spatio-temporal chaos', because their emergence results from collective action of turbulent 
fluctuations, as explained below. This phenomenon should be considered together with 
collective motion in high-dimensional chaos, which has been discovered in spatially extended 
systems [P|,pi][ and globally coupled systems [|r^,|r^. We believe that the collective patterns 



we consider presently are closely related to these types of collective motion. 

In this paper, we would like to focus on the universal aspects of the phenomena in 
question, not system dependent details. We first present a simple mathematical model 
which exhibits the emergence of collective patterns out of spatio-temporal chaos. We then 
explore how a transition occurs from a state of homogeneous spatio-temporal chaos to one 
displaying spatial order by analyzing this model. 

Equations of motion in fluid systems are written as 

dtq^'^ + = 0, (1) 

where q = {q^^\ ■ ■ ■ , g*-^^-*) is an M-component field defined in a D-dimensional space, and 
the current J = {J^^\ ■ ■ ■ , J^^^) is a functional of the field q, (e.g. q is given as g = {v, T) in 
thermo convective systems.) Eq.(|ip is often referred to as a continuity equation associated 
with some conservation law. We note, however, that Q^^^ = J dx^q^^^ is not a conserved 
quantity in a finite system where the surface integral of the current does not vanish. Nev- 
ertheless, when we focus on a bulk region in a turbulent system whose fluctuations develop 
only to a length scale much smaller than the system size, we can reg ard Q(^) nearly con- 
served quantity. This property plays an essential role in the formation of collective patterns, 
because large scale dynamics of conserved quantities determine the long-time behavior of 
fluctuations. Further, due to the fact that the wavelength of collective roll patterns is of 
the same order as the role height in convective systems, we conjecture that the formation 
of collective roll patterns is related to a breakdown of conservation laws. Based on these 
considerations, we propose the following one-dimensional model equation in order to study 
collective patterns: 



5tg^'^ + 9,/') = e'F«, (2) 

where F = {F^^\ ■ ■ ■ , F^^^) is a functional of the field q, and e is a small parameter which 
measures the extent of the breakdown of the conservation properties. The system size L 
is assumed to be much larger than the correlation length of fluctuations in spatio-temporal 
chaos, and a periodic boundary condition is assumed so as to maintain spatial translational 
symmetry. 

The simplest equation exhibiting spatio-temporal chaos with a conserved quantity is the 
Kuramoto-Sivashinsky (KS) equation, 

dtv + d^iv^ + d^v + dlv) = 0, (3) 

which has been derived in several contexts, including chemical waves, front propagation, 
etc. [0]. Statistical properties of turbulent solutions of the KS equation are homogeneous 
under periodic boundary conditions and have been investigated extensively Although 
it was found that the homogeneity of the statistical state is broken under a rigid boundary 
condition |T5|, we do not consider such an external symmetry breaking. Keeping in mind 
the guiding principle of modeling mentioned above, we studied a modified KS equation 
supplemented with a small term which breaks a conservation law. However, as far as we 
have checked, statistical properties seem to remain homogeneous if F is restricted to a local 
functional of v. Proceeding to the next step, we considered a two component system (i.e. 
M = 2) in which a second field is coupled to the KS equation, and found that the dynamical 
behavior is drastically changed. The model equation we discuss in this paper takes the form 

dtv + dx{v'^ + dxV + d^v) = —bv — du, 

dtu — Dd^u = au + cv — gu^, (4) 

where a,b, c, d and g are assumed to be small positive parameters proportional to e^. These 
values are chosen so that the trivial solution u = v = is stable in the absence of the term 
d^v, which produces spatio-temporal chaos. The condition for this stability is found to be 
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b > a and —ab + cd > 0, and we therefore assume throughout this paper that a = e^, b = 2e^, 
c = e"^, d = 3e^, and g = e^. 

Now, we present results of our numerical simulations of Eq.(^, where we have used a 
simple discretization scheme with a spatial mesh size 6x = 1.0 and time step size 6t = 0.1. 
This crude numerical scheme is sufficient for our purposes. (We note that our numerical 
scheme preserves a discrete spatial translational symmetry.) The values of e and L were 
varied between 0.1 and 0.025 and between 256 and 2048, respectively, because we are inter- 
ested in the system behavior under the asymptotic condition (e, L) (0, oo). The diffusion 
constant D is treated as a control parameter. In Fig. |l], spatio-temporal patterns are dis- 
played for different values of D, with (e, L) = (0.1, 256). In panel (a), small scale fluctuations 
are seen as usual in spatio-temporal chaos, while in panel (b), we can find that "a coarse- 
grained periodic pattern" is superimposed on the small scale fiuctuations. This suggests 
that our model can describe a transition to spatio-temporal chaos with a collective pattern. 
We confirmed quantitatively the existence of the average pattern by measuring the quantity 
< v{x) >, where we have introduced the abbreviated notation 

1 /-To+T 

<f>=7f; dtfit). (5) 

Here, Tq is some initial, and T is the averaging time. In order to check whether or not the 
field < v{x) > has a spatial variation in the limit T — oo, we measured the quantity 

(6) 



< v{x) > — — I dx < v{x) > 



= — f dx , . ^„ ^ , 
ij J ij 

where the spatio-temporal chaos is considered to be homogeneous when the scaling relation 
0" ~ 1/ \/T is obtained for sufficiently large T. Figure || shows the T dependence of a for 
systems in which (e, L) = (0.1, 1024). a decays to as a ~ 1/v^ when D = 1.5, while it 
seems to keep a constant value when D = 0.5. We thus expect that there is a critical value 
Dc between 0.5 and 1.5. The D dependence of a for a sufficiently large T [T = 12800), the 
phase diagram for our system, is shown in Fig. ^. This graph suggests that Dc is near 0.8. 

We now attempt to explain how the collective patterns emerge out of spatio-temporal 
chaos. First, let us derive a self-consistent condition under which a long time averaged 



quantity < v{x) > forms a periodic pattern. For simplicity, we focus on the asymptotic case 
that (e, L) (0, oo) and suppose the e dependences of < m > and < v > can be expressed 
as < M >= U{ex), and < v >= V{ex). These scahng forms are consistent with numerical 
results. Then, substituting u =< u > +u' and v =< v > +v' into Eq.(^), averaging the 
equation, and discarding o(e^) terms in the equation, we obtain 

d^{V^+ < v'^ > +d^V) = -bV - dV, 

-DdlU = aU + cV - gU\ (7) 

where we have assumed < {u— < u >)^ >= o(e) with = 2, 3. Now, in order to obtain a 
closed form of the equation for U and V , we employ the phenomeno logical relation 

< v'"^ >= a - vd^ <v >, (8) 

which has been used to analyze turbulent shear flow 0. Here, a and v are assumed to 
be constants, and v is often referred to as "turbulent viscosity". We can see from Eqs.(^ 
and that v expresses how turbulent fluctuations contribute a dissipation effect on the 
coarse-grained pattern. Although the validity of Eq.(^) for our equation is not proved 
mathematically, as shown in Fig. ^, our numerical results for < f'^ > and < v > are 
consistent with this relation. Further, under a mode truncation approximation, we write 
< V >= Rsm{kx + (p) and < u >= S sm{kx + 0), and obtain the equation for R and 5" 

{D,k'^ + b)R + dS = 0, 

-cR + {Dk^ - a)S = -gS\ (9) 

where k is of order e, in accordance with the above stated assumption, and D^, = v — 
1. Equation (^) always has the trivial solution i? = 5* = 0, which represents a spatially 
homogeneous statistical state, while a non-trivial solution exists only when the inequality 

D < '^-^D^ (10) 

is satisfied. This result is consistent with our phase diagram shown in Fig.3. In fact, by fitting 
our numerical data to Eq.(|]), we found v ~ 5.5 in the case that (D, e, L) = (0.2, 0.1, 1024). 
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Using this value of v for all we obtain ~ 0.6. The discrepancy between this value and 



The self-consistent approach developed above can tell us nothing about the stability 
of the obtained solutions. In order to consider the stability problem, we need to develop a 
theory describing the slow evolution of coarse-grained patterns. A derivation of the evolution 
equation from our model equation is beyond the scope of this paper, but we present the 
results we obtained using the scaling relation such that f (x, t) = V{ex, eH) + t). We 
found that the homogeneous statistical state becomes unstable when the inequality (p!0D is 
satisfied, and in this case, a coarse-grained periodic pattern with a wavenumber grows 
with the maximum growth rate, where 



That is, a non-trivial persistent pattern appears when the homogeneous statistical state 
becomes unstable. This statement is consistent with our numerical experiments. Here, we 
would like to stress that D^. and are determined by D^,, which is not a parameter of 
our model but associated with correlations of fluctuations. Therefore, the existence of the 
averaged pattern cannot be explained by developing a linear stability analysis of the trivial 
state u = V = 

As argued above, the collective action of fluctuations causes an instability leading to the 
formation of persistent patterns. In this sense, we can say that fluctuations are the origin 
of a spatial order. On the other hand, one may suspect that fluctuations act as a turbulent 
noise, inducing a random drift of our coarse-grained pattern, like Brownian motion, and as 
a result the averaged pattern vanishes in the limit that the averaging time T goes to infinity. 
In order to check this, we calculated a time series of the complex amplitude of the Fourier 
mode associated with the collective pattern, {A{t) exp(i0(t))}. Spectra of the amplitude A 
and phase 0, which are denoted by Sa{oj) and S^iu) respectively, are shown in Fig. ^. The 
spectrum S^ioj) goes as u;"^ at low frequency {uj < ujc). This implies the existence of a phase 





(11) 
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diffusion of the coarse-grained pattern and therefore that the averaged pattern will vanish 
when the averaging time T is chosen much larger than 2tt/uJc- We should note, however, 
that the crossover frequency Uc is much smaller than the inverse of the correlation time of 
chaotic fluctuations, and from our preliminary numerical results we believe that tUc — > for 
e — ^> 0, with tL fixed. This implys that if the limit T — oo is considered after taking e — > 
with fixed tL > 0(1), persistent patterns actually exist. Figure 5 shows further that the 
amplitude R has a long-range correlation. At low frequency {u < Uc), the form of Sr^u) 
seems to be a Lorenzian type, proportional to (1 + (cur)^)"^. Here, the correlation time 
of the amplitude fluctuation, r, which may depend on (e, L), is much larger than 2tt/uJc- 
Although we believe that tuJc — oo for e ^ with fixed eL, it is necessary to study the 
(e, L) dependence of the system behavior further. 

In conclusion, we have shown that our model equation can exhibit the emergence of 
collective patterns out of spatio-temporal chaos. We believe that such a phenomenon will 
be observed in other mathematical models which have nearly conserved quantities. By 
introducing a small parameter e which measures the extent of the breakdown of conservation 
properties, we have developed a self-consistent analysis which describes a transition to the 
spatially ordered phase. Both fluctuations of the amplitude and phase of the Fourier mode 
associated with the collective pattern have long range correlations at a time scale larger 
than 2ti/uJci while over a time scale 211 /uc-, which goes to infinity for e — > 0, averaged 
patterns are properly defined. In our theoretical analysis, the relation between the averaged 
pattern and the correlation of fluctuations around it plays a crucial role. Although the 
relation (P) is a familiar one at a phenomenological level, there may be various forms which 
relate the statistical property of fluctuations with an averaged quantity. We expect that as 
yet unknown relations will be discovered in experimental systems and other mathematical 
models. Also, to determine the correspondence between our model and experimental systems 
may be an important subject. Such studies should contribute greatly to the understanding 
of the nature of spatio-temporal chaos. 
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FIGURES 

FIG. 1. Spatio-temporal patterns of the v field of our model equation with (e, L) = (0.1, 256). 

We have, (a) D = 0.5 and (b) D = 1.5. The time goes from t = 10000 to t = 22800. 

FIG. 2. a versus T for L» = 0.5 (solid line) and D = 1.5 (dashed line) (e,L) = (0.1,1024). 
The transient time was chosen as Tq = 10000. 

FIG. 3. a versus D for Tq = 10000 and T = 12800. (e, L) is as in Fig. 2. Initial conditions 
are the same for all values of D. 

FIG. 4. Snapshot of the v field at t = 505000 (dotted line), the averaged pattern (solid line), 
and the dispersion at each point (dotted line), for the system with {D, e, L) = (0.2, 0.1, 1024). Each 
field in the region < x < 256 is displayed. The statistical procedure was carried out for a sample 
{v{x,ty,5000 < t < 505000}. 

FIG. 5. Spectra for a time series of the amplitude (solid line) and phase (dotted line) of a 
Fourier mode exp{ikx) associated with a collective pattern with the same parameter values as in 
Fig. 4. The value of k is given as the position of the peak of the spectrum of the averaged pattern. 
These graphs were obtained by averaging 8 samples of the spectra using 8192 data points every 
128 time units. 
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Fig. 1a S. Sasa 




Fig. 1b S. Sasa 




Fig.2 S. Sasa 




Fig. 3 S.Sasa 




Fig.4 S.Sasa 




Fig. 5 S.Sasa 



